subroutine get_estop(e,se)
  use param
  use global
  implicit none
  real(dbl),intent(in) :: e
  real(dbl),intent(out) :: se
  se=k4*e**0.5
  return
end subroutine get_estop
